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Abstract. 

Four families of generalizations of trigonometric functions, denoted cos^ ^j(a;, y) and sin^ (a;, j/), 
were recently introduced. In the paper the functions are transformed into four families of orthogonal 
polynomials depending on two variables. Recurrence relations for construction of the polynomials are 
presented. Orthogonality relations of the four families of polynomials are found together with the appro- 
priate weight fuctions. Tables of the lowest degree polynomials are shown. Numerous trigonometric-like 
identities are found. 

Two of the four families of functions are identified as the functions encountered in the Weyl character 
formula for the finite dimensional irreducible representations of the compact Lie group 5*^(4). The other 
two families of functions seem to play no role in Lie theory so far in spite of their analogous 'good' 
properties. 
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^ 1. Introduction 

^ ^ ^ Four families of special functions and four related families of orthogonal polynomials in two real 
variables are studied in the paper. They are the symmetric and antisymmetric generalizations of cosine 
and sine functions denoted here as cos^ and sin^. The simplicity of the functions ([T]) combined with 
their continuous and discrete orthogonality should make them particularly amenable to application. 
Orthogonal polynomials were studied intensively in recent years using various approaches frequently 
10 inspired by orthogonal polynomials in one variable. Many such polynomials are described in the liter- 
ature [H [21 |3l H] . In the absence of additional constraints, truly higher dimensional generalizations are 
^— I hidden in the vast number of possibilities of defining orthogonal polynomials of more than one variable 
121 El m [7] . In 2D [1] , such constraints were provided by requiring that polynomials of two variables be 
simultaneous eigenvectors of two differential operators. 
^ The idea of our construction of the polynomials in its simplest realization is rater transparent: Con- 
• ^ sider the characters of the finite dimensional irreducible representations of a compact simple Lie group 
^ of rank n. Declare the characters Xuii , ■ ■ ■ , Xujn of the fundamental representations to be the polynomial 



^ variables. Any other character is then a polynomial in the lowest n of them. It is constructed recursively 
by a judicious choice of products of the lower characters with the n basic ones and their decomposition 
into the sum of irreducible characters. Orthogonality and other properties of such polynomials follow 
from the corresponding properties of characters. 

The idea of our construction remains valid when the characters are replaced by their constituents. 
That is the Weyl group invariant orbit functions (C-functions here), or the Weyl group skew-invariant 
orbit functions (5-functions here) that appear in the numerator of the Weyl character formula, or the 
'hybrid' functions C~ and defined in this paper. 

Our starting point are the special functions of n real variables [8], which we specialize to two variables, 
and from those functions we construct the 2-variable polynomials. In order to do that, we substitute 



the variables (16) into the functions of (14). The functions thus become polynomials in X and Y. The 
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recurrence relations of subsections 4.2- 4.5 allow one to find the polynomials to any desired degree. 
Several lowest degree polynomials for each of the four families are shown in Appendix 2. 

The immediate goal of our paper is to describe properties of such functions and their polynomials 
in details impossible in the n-variable generality. In this way we obtain only a limited number of the 
polynomials of 2 variables found in the literature, but our approach has other advantages besides the 
possibility of generalization to any number of variables obvious from ^ . Let us now point out the main 
advantages. 

The most important conclusion of the paper is the link between the special functions, generalizations 
of the symmetric and antisymmetric sine and cosine, and their polynomials to the compact simple Lie 
group of type C2 (also called 0(5) or 5^(4)). Relating the polynomials to the compact simple Lie 
groups, rather than to the root systems of such groups [1], provides a discretization of the polynomials 
through the discretization of the special functions [9]. Discrete orthogonality of the functions carries 
over to discrete orthogonality of the polynomials on lattices of any density and of symmetry dictated 
by the underlying Lie group. See examples in \10\ llll [T2] . Moreover a generalization to the simple Lie 
group of any rank > 2 is straightforward and leads to polynomials of any number of variables |10| . 

The most curious result of the paper is the existence of two new families of orthogonal special 
functions, called here C~ and 5"*", besides the families of C and S types which are known from the 
Weyl character formula [13j. Analogous new families of special functions are to be studied elsewhere 
for every compact simple Lie group with two root lengths, the number of variables being the rank of 
the Lie group. The new families of functions do not exist for the simple Lie groups with roots of equal 
length, i.e. types An, Dn, Eq, £'7, and Eg. 

Previous results |14|, [T5] about the C2-functions are extended in this paper in several ways: 

- One-to-one correspondence is shown between functions of the families cos"^, sin~ and the C- and S- 
functions appearing in the Weyl formula for characters of finite dimensional irreducible representations 
of the compact simple Lie algebra C2 of the Lie groups 5p(4) and 0(5). The correspondence is a result 
of a straightforward linear substitution of variables The fundamental weight basis and its dual in 
the case of C- and /S- functions are transformed into the orthogonal basis of cos"^ and sin~ functions. 

- The same substitution of variables reveals that the functions cos~ and sin"*" correspond to two families 
of functions, denoted here as C~ and S^, whose role in Lie theory has apparently not been noted before. 
The functions of the new families have all the 'good' properties of the C- and S-functions, which are 
called here the C~^- and ^"-functions. The notable distinction between the four families of functions is 
their behaviour at the boundary dF of the isosceles right triangle F, the domain of their orthogonality 
in the real Euclidean space M^. 

- A method of constructing polynomials of two variables from the special functions is developed and 
applied to finding some lowest degree polynomials for each family. It can be understood as a non-linear 



substitution of variables linking F C M and its image 5^ C M after the substitution of variables (14) 



Both F and ^ are shown in Fig. 6.1 



- Continuous orthogonality of the polynomials, when integrated over J, is demonstrated. An appropriate 
weight fuction is determined for each famil}|^ It is a consequence of the orthogonality of the irreducible 
characters of any compact simple Lie group when they are integrated over the maximal torus, and of 
the symmetries of our functions under the action of the Weyl group. 

- Numerous trigonometric-like identities are found involving generalized trigonometric functions of the 
four families. Many of them are shown in Appendix 1. 

The history of the generalization of trigonometric functions considered here can be traced in the 
three recent papers [HI [HI IS]- The C-functions of the simple Lie algebra C2 were considered in [TUfTB]. 
and the 5- functions of C2 are described in [ 12| [T7] . The polynomial representation of the functions ([T]) 
has not been considered elsewhere. In general, the study of the orthogonal polynomials of C-and S- 
functions originated in [18]. Their exploitation in polynomial approximations [19j is another possibility. 
Simultaneous zero values of polynomials of certain degrees have to be identified and the corresponding 



"'^The discrete orthogonality of the polynomials, when summed up over a finite set of points, is described elsewhere for 
each family and for several low rank simple Lie groups. 
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cubature formulas used [lOl EQl [5]. Computational aspects of such applications remain unexplored 
besides the few examples in [HI [I5j . 

In Section 2, the functions cos^^^(x,y) and sinp^^^(x,y) are introduced and their congruence class 
is defined. The structure of the decomposition of their products is shown in Table [l} Changing from 
the orthogonal basis to bases inherent to the C2 symmetry in Section 3 shows that the functions ([T]) 
are in fact, up to a normalization, the VK-invariant functions of the simple Lie group C2. Two of 
the functions are well known in Lie theory, the other two are new. The polynomials are described in 



Section 4. The functions (14) to be represented as polynomials in the same variables are considered 
first. Curiously, the structure of three of the four families of functions resembles the Weyl character 
formula. In fact using the variables of the previous Section, it is easy to verify that the function sin~ is 
the character function of C2. The variables X and Y are then introduced, the lowest few polynomials 
of cosp^ (x, y) are constructed, and the full set of recurrence relations is written for each family. The 
functions studied here are shown in Section 5 to be eigenfunctions of the Laplace operator and of another 
algebraically independent 4-th degree operator. In orthogonal coordinates this fact is visible from ([T]). 
In the polynomial variables both operators are considerably more complicated. Section 6 examines the 
continuous orthogonality of the polynomials, known for the functions ([T]), and rewritten here in the 
basis {X, y}. Unexpectedly, the kernel-Jacobian of the orthogonality integrals known from Lie theory 
factorizes if it 'knows' of the existence of the two new families of functions. Section 7 contains the 
generic identities verified by the functions ([T]). Special identities are shown in Appendix 1. The last 
Section contains some concluding remarks. Appendix 2 consists of eight tables of explicit low degree 
polynomials, two tables per family. A comparison of the lowest polynomials between different families 



of (14) can easily be made there. 

2. Notations, definitions, and some properties 

We first recall the definitions of functions symmetric and antisymmetric in variables x and y, intro- 
ducing simplified notations, 

(x, y) = cos(7rAx) cos(7r;uy) -|- cos(7r^x) cos(7rAy) (symmetric cosine) 
(x,y) = cos(7rAx) cos(7r;uy) — cos(7r^x) cos(7rAy) (antisymmetric cosine) 

(1) 





:= cos,\ s 




"°^(A,M) 


+ 

(A,/.) 


^i^fA.M) 




'■= ^i^rA,M) 



(x, y) = sin(7rAx) sm(TTfiy) + sin(7r/ix) sin(7rAy) (symmetric sine) 

(x, y) = sin(7rAx) sin(7r^y) — sin(7r/ix) sin(7rAy) (antisymmetric sine) 
Here x, ?/ G M and A, G Z, such that 

A>/i>0 for c+ _^), X>fi>0 for c^^ 

A>/i>0 for A>^>0 for s'^^^^ 

Note that c~ and s" are 2x2 determinants, while c"*" and are 2x2 permanents [21j. Assuming that 
A and are integers, it is advantageous to split the functions within each family into two congruence 
classes according to the value of the congruence number # defined as 

#(c±_^)(x,y)) = #(5±_^)(x,y)) = A-^ mod 2 . (3) 

During multiplication of the functions their congruence numbers add up. 

Numerous symmetry properties of the polynomials Cp^ and are shown in [TJ] and [IS] re- 

spectively. The functions ([T]) are clearly defined for all real values of x and y and have continuous 
derivatives of all degrees. Due to their symmetry properties, we are interested in their values only in 
the fundamental region F cM?, which is the right isosceles triangle on Fig. 6.1 with vertices 



F:={(0,0), (1,0), (1,1)}. (4) 

Of importance is the distinct behaviour of functions ([T]) at the boundary dF of F. Indeed, normal 
derivatives of c"*" functions vanish at dF, functions s~ are equal to zero at OF, while c~ and behave 
differently at long and short sides of the triangle F. 

Products of pairs of the functions ([T]) decompose into their sums. 10 different pairs of functions can 
be multiplied. The decompositions have a rigid structure, which affects which functions may appear. To 
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demonstrate, we reproduce Table 1 of [15] as Table 1 here. Actual decompositions of specific products 
are found in Table 2 of |15j . 



product 


sin+ sin+ 


sin+ sin 


sin sin 


sin+ cos"*" 


sin"*" cos 


sin COS+ 


sin cos 


cos"*" cos"*" 


COS+ cos 


COS COS 


terms 


cos"*" 


COS" 


cos^ 


sin^ 


sin~ 


sin~ 


sin^ 


cos^ 


COS" 


cos"*" 



Table 1. Structure of the decomposition of the 10 types of products. The second row 
shows which functions appear in all the terms of a decomposition. 



3. Orbit functions of C2 and cos='= and sin='= 

The C- and S'-functions defined by the summation of exponential functions over the Weyl group 
W{C2) were recently studied as special functions with many practically useful properties. In Lie theory 
they are known from the Weyl formula for the characters of irreducible representations of the compact 
Lie group C2. We show that C-functions coincide (up to a scaling factor) with c"'"-functions defined 
here, and that S'-functions similarly coincide with ^"-functions. Moreover, we find analogs of c~ and 
that were not noted in Lie theory before. We describe their properties in parallel with the more 
familiar C- and S'-functions. See [22j for a study of general simple Lie groups not only of type C2. 

An orbit of W{C2) consists of all the distinct points (weights) obtained from the dominant point v_ 
by repeated application of reflections Tq,^ and in mirrors orthogonal to the simple roots ai, a2 of 
C2, according to the formula 

2(a,-,v) 

ra.V = V- rUj . 

{aj,aj) 

Assume that the qi is the shorter of the two simple roots of C2, i.e. 

(«i,ai) = l, (02,02) = 2, (ai,a2) = -l- 

In this paper we set v = {vi,V2) = viUJi + ^2^2, (fi,f2 S Z-^), and 9 = (^1,^2) = ^I'^i -|- O2UJ2, 
(^1,^2 £ Recall that the scalar product of two vectors, one given in the cj-basis and the other in 
the dual tj-basis of C2, is calculated as follows |16| . 



2 

iV.,^) = ((^"1^1 + ^2^2), (6*1^1 + 92UJ2)) = ^ ViOjioJi^Lbj) = {Vl + V2)0l + (i-Ul + V2)02 ■ 

i,j = ^ 

A generic VF(C2)-orbit consists of 8 points/weights: 

±{vi,V2), ±{-vi,vi + V2) , ±(^1 -I- 2^2, -f2) , ±(t!i -I- 2^2, -fi - ^2) ; t'l,f2GN. 

The generic C- and S-functions are denoted in the sequel by C"*" and S~ respectively. Together with 
the new functions C~ and S"*", they are defined and identified with the symmetric and antisymmetric 
trigonometric functions of 2 variables as follows. 

(5) 
(6) 
(7) 
(8) 





= E 




= 4cJ„^)(a;,y), 




= E 

wew(C2) 


■ _-^yi{w)+l2{w) ^2ni(wv,9 






= E 

wew(C2) 








= E 

w£WiC2) 







Here hiw) is the minimal number of times the reflection r^^. is required to generate w. 



The equalities of the functions in ([5]) - ([s]) are demonstrated using the fohowing change of uj- and 
w-basis to the orthogonal basis {61,62}, where (ei,6i) = (62,62) = 5, and the corresponding variables 
vi, V2 are substituted for A, /i, and 61, 62 for x, y, 



wi = 6i, wi = 26i, fi = A-/i, 0i = ^{x-y), 

W2 = 61 + 62 , W2 = 61 + 62 , V2 = 02 = y, 

In particular, #v = #{vi,V2) = vi mod 2. 

Using our new symbols for C-functions and for S'-functions, we have explicitly, 

C+^,,,)(^i,^2) := C(,,,,,)(ei,^2) = 2cos{27r{{vi,V2),e)) + ^cos{27T{{-vi,vi+V2),e)) 
+ 2cos(27r((wi + 2t;2,-U2),^)) + 2cos(27r((wi + 2r;2,-ui - ^2),^)) 
= 2 (cos(7r((i;i + t>2)26'i + {vi + 2i;2)6'2)) + cos(7r(w226'i + {vi + 2t;2)6'2)) 
+ cos(7r((wi + 7;2)26'i + W16I2)) + cos(7r(t;22(9i + {-vi)e2))) ■ 



(9) 



(10) 



fill 



5(,,,,,)(^i,^2) := 5(,,,.2)(^i,^2) = 2cos(27r((t;i,t;2),^)) - 2 cos(27r((-t;i, t;i + i;2), 0) 

- 2 cos(27r((t;i + 2^2, -^2), ^)) + 2 cos(27r((?;i + 2^2, -vi - V2),0)) 
= 2 (cos(7r((t;i + r;2)26'i + (vi + 2v2)62)) - cos(7r(?;22(9i + (vi + 2^2)612)) 

- cos(7r((t;i + V2)2ei + ^16*2)) + cos(7r(t>226'i + (-fi)6'2))) • 

Clearly the C"*" and S~ differ by the sign of some terms. More precisely, the sign is negative each time 
an odd number of reflections r^^, r^j is applied to v = (^1,^2)- A practically important difference 
between functions of the two families is their behaviour at the boundary dF of F. The function is 
symmetric with respect to reflections across the three sides of F, while S~ is antisymmetric. Both are 
continuous differentiable, therefore normal derivatives of at dF are zero, and has zero value at 
OF. 



The fundamental role of the functions (10) and (11) in representation theory is asserted in that they 



comprise the Weyl character formula for the Lie group C2. 

The substitutions ^ used in the functions C~ and result respectively in cos~ and sin+, which 
do not play any noted role in Lie theory, as far as we know. 

'^(~.i,.2)(^i'^2) := C(,,,,2)(^i>^2) = 2cos(27r((t;i,t;2),^)) - 2 cos(27r((-?;i, + t;2), 0) 

+ 2cos(27r((vi + 2t;2,-i;2),^)) -2cos(27r((wi + 2T;2,-7;i -^72),^)) 

= 2 (cos(7r((vi + ^2)2^1 + {vi + 2v2)92)) - cos(7r(w22^i + {vi + 2^2)^2)) 

+ cos(7r((wi + 7;2)26'i + W16I2)) - cos(7r(t;22(9i + (-vi)(92))) . 

Negative signs of terms in (|12[) occur whenever an odd number of reflections r^j is applied tov = (^1,^2)- 



5j,^,,,)(ei,^2) := 5(,,,,,)(ei,^2) = 2cos(27r((t;i,t;2),^)) + 2cos(27r((-t;i,t;i +^2),0) 

- 2 cos(27r((t;i + 2v2, -V2),0)) - 2 cos(27r((t;i + 2^2, -vi - V2),0)) 

= 2 (cos(7r((t;i + ^2)2^1 + {vi + 2v2)92)) + cos(^(t;220i + {vi + 2^72)^2)) 

- (cos(7r((t>i + V2)2ei + vi92)) - cos(7r(t;226'i + (-wi)6'2))) . 

Negative signs of terms in (|13[) occur whenever an odd number of reflections is applied tov = (vi,V2). 



4. Construction of polynomials of cos='= and sin='= 

A polynomial is irreducible if it belongs to one of the four families ([T]) and if it carries specific values 
of the integers A and /x, subject to the inequalities ^ applicable to its family. 
Our aim in this section is to represent the following functions, 

cos7\ Ax,y) sint, ^(x,y) sin7, Jx,y) 

COS+ {x,y), (^'^)^ ^ ^ __(M^)\^yj_ ^ _JM^y;_^yj_ ^ A, //GZ-°, (14) 

' ' cos^^_Q^(x,y) ' sin+^-^^(a;,y) ' sin^2,i)(2;' y) ' 
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as polynomials in the variables X and 1" of (16). In addition, one must require that A and /i satisfy the 
inequalities prescribed in ^ for each type of function. Furthermore, if we multiply functions and 
s~ by certain constant gx^ given by 



3\n 



1 A = ^ = 

2 A = ^ > 
2 A > ^ = 
4 otherwise 



(15) 



we obtain orbit functions defined by the summation of exponential functions over the orbit of the 
Weyl group instead of the summation over the Weyl group |16| [T7] . In [18 1 , there are constructed 
polynomials of such orbit functions, such polynomials are monomials. Thus, we study here the functions 



( PI multiphed by c/a,/.- 

By direct computation, one may verify that the three functions in the denominators of ( 14 ) have no 



zero values inside the fundamental region F of C2. Moreover, whenever one of the denominator functions 
turns out to be zero at the boundary dF, the family of corresponding numerator functions is also zero 
there. More precisely, the functions c~ are zero at the long side of the triangle F, and the functions s"*" 
are zero at the short sides of F. Such behaviour is known in general only for Sg, where g is the half sum 
of positive roots of any compact simple G. In our case, —4 sinT^ u) = ^dx, y) = S*)^! i)(x, y). 



Our main tool in constructing the polynomials is the decomposition of products of pairs of functions 
([T|, besides the choice of variables. Since all products were decomposed in [1^, we refrain here from 
describing methods of finding the decompositions. 

Following the general idea of [18j, the lowest functions multiplied by 2 serve as the polynomial 



variables X and Y for all four cases (14): 



X := 2c^ y) = 2(cos(7rx) + cos(7ry)) , y := 2c^ ^^(x, y) = 4 cos(7rx) cos(7ry) , x,yGM. (16) 



(1,1)^ 



Clearly we also have n~i(^' ~ ^- addition, we need the three variables appearing as denominator 

(17) 



terms in (14): 



"(0,0) 
U := 2c 



(1,0) 



[x,y) , 



y :=2s+,)(x,y) 



W := 4s(2^-^^(x,y) . 
#W = 1, 



4.1. The polynomials cos"*". 

According to Table [l| polynomials of type arise in decompositions of products of four different 
types, namely cos^cos"^, cos~cos~, sin"*" sin"*" , and sin~sin~. Decompositions of products X'^ and 
yield quadratic polynomials, both with the congruence number # = 0. The only polynomial quadratic 
in X and Y with 7^ = 1 is XY. Calculating the products of two variables, we obtain 

2c+ 0) =X^ -2Y -A, 2c+ 2) =1^^-2X^ + 4^ + 4 

Taking ever higher functions multiplying them by either X or Y, and decomposing the products, 

we obtain higher and higher degree polynomials. Continuing with polynomial of # = 0, we obtain the 
next lowest polynomials, 

40^3 ^ = X'^Y - 2y2 _ QY 

2c+ 3) = Y^ - 3X^Y + 6Y^ + 9Y 

2c+ ^^ = X^- AX^Y + 2y2 _ + sy + 4 ^^g^ 
4cJ 2) = ^^^^ - 2^*^ - 2y^ + BX^y - 12y2 + lOX^ - 20Y - 8 
2c+ = y4 _ 4x2y2 + 2X'^ + sy^ - sx^y + 20y^ - sx^ + ley + 4 

See Table [| for more c"*" polynomials with # = 0. 
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Decompositions of the products with ^ = 1 are obtained in a similar way: 



4ci, ,^ 
(2,1) 


= xy - 


- 2X 


^^(3,0) 




3XY - 3X 


4r+ 
*C(3,2) 


= 


- 2X^ + 3XY + 6X 


*C(4,1) 




- 3Xy2 _ 4XY + 2X 


*''{4,3) 


= xy'^ 


- 3X^Y + 5Xy2 + 2X 



(19) 



6X 



Note that there is an easy useful check on any recurrence relation or decomposition of a product. Indeed, 
specializing {x,y) = (0,0), we get 

X(0,0) = y(0,0) = 4, (20) 



(0,0) 



which turns any recurrence relation into an equality of integers. 

See Table [2] and Table |3] for more polynomials with # = and 1 respectively. 



4.2. Recurrence relations for symmetric cosine functions c 



(A,/x)- 



Using the variables (16) as before, and decomposing appropriate products of c"*" functions, we have 

2c1 



A > 2, /u = : 


2c+ 

(A,0) 


= 4c 


A - 1 > ^ > : 


(A,m) 


= 8c 


A - 1 = ^ > : 


4r+ 


= 4c 


A = ^ > 2 : 


2c+ 

^^(A,A) 


= 4c 


Again the check (20) is valid 


in (21). 





-(A-1,/.)^(1,0) 
-(A-1,A-1)''(1,0) 
(A-l,A-l)C(l,l) 



4c 



'•(A-2,0) 
+ 

(A-2,m) 

4c' 



4r+ 

(A-1,1) 



4c 



(A-l,/x+l) 



4c: 



(A-l,/.-l) 



(21) 



4c 



(A-1,A-2) 

+ 

(A,A-2) 



2c 



(A-2,A-2) 



Starting from the lowest admissible values of A and we recursively build c"*" polynomials up to any 



desired degree in terms of variables X and Y using (21). Some insight into the structure of recurrence 
relations is gained when they are rewritten as 3-term relations with matrix coefficients [6l |7] . 

Introduce the polynomials renormalized by the factor (15) := ^A^tCj^^-j. Form the vector Fx, 

Pa := (Pa,o(^, Y),Px,i{X, Y),..., Pa,a-i(^, Y), Pa,a(^, Y)f G M^+i'^ (22) 
and define on Pa the operation 

{X, Y)Fx := (XFI, YFlf G M^^+^.i (23) 

Assuming that the lowest polynomials with A < 3 have been calculated, see Table [2] and Table [3j the 
polynomials with A + 1 > 4 are found from the 3-term recurrence relation: 

Pa+1 = Da(^, :^)Pa + EaPa + FaPa-1 , (24) 
where the matrix coefficients are the following, 

/I 



V -1 



G 



i)A+2,2A+2 



/-I 




1 



V 

/O -1 
-2 -1 

-1 -1 



G 



-1 -1 

-1 -2 
0-10 

1 2 / 



G 



pA+2,A+l 
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4.3. Recurrence relations for antisymmetric cosine functions C^Xfj.)' 

Polynomials of type c~ arise in decompositions of two kinds of products ([T]), namely cos^cos", and 
sin^sin". Recurrence relations for the function are described in terms of the variables X, Y, and 



U, where hU = c,, „n is the lowest of the c -functions. 



2^ "-(1,0) 

There are two second degree polynomials 



2c 



(2,0) 



ux , 



4.C 



(2,1) 



U{Y + 2) 



Decompositions of appropriate products yield the recurrence relations for higher degree polynomials. 



X>2, fj, = 
X-2> fi>0 
X>2, X-l = fi>0 
X-2 = fi>0 



^"^(A.O) ~ ^'^(A-1,0)'^(1,0) ^'^(A-2,0) 

^C(A,M) = ^^(A-l,^)C(l,0) ~ ^^(A-2,m) 

^•^(A.A-l) ^ ^'^(A-1,A-2)'^(1,1) ~ 

^'^(A,A-2) ^ ^'^(A-1,A-2)'^(1,0) ~ ^'^ 



Ac 



(A-1,1) 



4c 



(A-l,At+l) 



4c 



(A-l,At-l) 



(A,A-3) ^'^(A-l,A-2) 
(A-l,A-3) 



4c: 



(25) 



(A-2,A-3) 



Using (25) and starting from the lowest admissible values of A and ^u, we build the c polynomials 



up to any desired degree in X, Y, and U, where U occurs in the first power only. Instead it is more 
interesting to construct the polynomials c^^^/C/ in order to compare them with the other polynomials 
of the functions ( |14[ ), all in X, y variables. 

Recurrence relations can be written in their 3-term form for the polynomials and solved provided the 
polynomials ^^/C/ for A < 4 are known. See Table jij and Table 5 



Introducing the vector ¥x with polynomial components, 

Pa := {Px,o{X, Y),Px,i{X, Y),..., Px,x-2(X, Y), Px,x-i{X, Y)f G M^-^ 
and the transformation 

{X,Y)Fx := {XFI,YFI)'^ £ M^A,! 
we get the recurrence relations for A + 1 > 5 

Pa+1 = Ox{X, Y)Fx + ExFx + F^Pa-i 



(26) 
(27) 
(28) 
(29) 



where 



/I 



/-I 



1 

-1 



pA+l,2A 



V 



-1 



0/ 



pA+l,A-l 



/O -1 

-2 -1 
-1 



-1 

-1 



-1 



-1 -1 
0-10 
1 2 / 



pA+l,A 



4.4. Recurrence relations for symmetric sine functions s 



Polynomials of type arise in decompositions of two kinds of products ([T]), namely cos^sin^, and 
COS" sin~. Recurrence relations are constructed for polynomials s^x i^) variables X, Y, and V, where 
= s^ii) is the lowest of the ^'''-functions. 

The quadratic polynomials being 4s^ = VX and 2s^ 2) = : the higher degree polynomials are 
found from the relations 



A > 2, 




= 1 : 


*^(A,l) 


= 


c+ - 

-1,1) (1,0) 


^*(A-2,1) 


4s+ 

(A-1,2) 


A-1 > 




> 1 : 




= K~ 


^i,m)^i,o) - 


4s+ — 


4s+ _ 4s+ 


A-1 = 




> 1 : 


^*(A,A-l) 


= 4«fA- 


-1,A-1)''(1,0) 


— 4s+ 


-2) 


A = 




> 2 : 


2s+ 
(A,A) 


= 


c+ 

-1,A-1)''(1,1) 


^'^(A,A-2) 


- 2s+ 

'^*(A-2,A-2) 



(30) 



As in the previous subsection, we are more interested in polynomials s^/F, which depend on X and Y 



only and therefore can be compared with polynomials of types (14). 
See Table joj and Table jrj of polynomials Sp^^^/y. 
Putting 



5'A,M'5(A,/i) 



/V 



and assuming that the polynomials are known for A < 3, see Table [6] and Table [7| we introduce the 
vector W\ as before, 



Pa := {Px,i{X, Y),Px,2{X, y), . . . , Pa,a-i(^, Y),Px,x{X, Y) f G 
together with the transformation 

{X,Y)Fx := {XFI,YFI)'^ £ M^^'^ , 
we have the reccurence relations for A + 1 > 4: 

Pa+1 = DxiX, Y)Fx + EaPa + FaPa-i • 

Here 

/-I 



i>A,l 



/I 

V 



1 

-1 



V 



pA+l,2A 



Fa 



(31) 
(32) 
(33) 



[,A+1,A-1 



/O -1 

-1 -1 



-1 



-1 



-1 -2 
0-10 
1 2 



pA+l,A 



4.5. Recurrence relations for antisymmetric sine functions s 



Polynomials of type s arise in decompositions of two kinds of products ([T]), namely cos^ sin , and 
COS" sin+. The variables in this case are X, Y, and W, where = is the lowest of s~-functions. 
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(3,1) 



WX and 4s 



(3,2) 



W{Y + 1) are then used together with 



4^(A,l) 


= 8^(A- 


r+ - 
-1,1) (1,0) 


4^(A-2,1) 


4S(A- 


-1,2) 


^^("a,m) 


= 8^(A- 




4^(A-2,M) " 


4^(A 




(A,A-2) 


= 


-l,A-2)^(l,0) 


~ 4'^(A-1,A- 


-3) 




(A,A-1) 


= 


-l,A-2)^(l,l) 


+ 4S^A-1,A- 


-2) ~ 


"^•^(A.A-S) ~ '^'^(A-2,A-3) 



(34) 



The two quadratic polynomials 4s 
the relations 

A > 3, /X = 1 
A-2 > /X > 1 
A-2 = ;U>1: 4s 
A>3, A-1 = ;U>1: 4s 

As in the previous subsection, we are interested in polynomials s~ /W which depend on X and Y only 
and therefore can be compared with polynomials of types (14). 
See Table jsj and Table ojof polynomials s^^^/VF. 

In the case s^ ^)/^' polynomials 

Px,^, = <?a,mS^a,^)/^ 

are suitably normalized irreducible characters of the simple Lie group C2. Introducing the vector 

\T ^ Tn,A-l,l 



Pa := (PA,l(^,^),^A,2(^,n,---,^A,A-2(X,n,^A,A-l(^,>^))' G 
and the transformation 



{X,Y)Fx 

we get the recurrence relations for A > 4: 



{xfI,yf 



G 



!)2A-2,1 



where 



Vi = Bx{X, Y)Fx + EaPa + FaPa-1 , 

/-I 



(35) 
(36) 
(37) 



/I 




\ 












G ]rA,2A-2 ^ 


Fa = 




1 






v 


-1 


V 










/o 


-1 








-1 


-1 










-1 -1 






Ea = 














-1 


-1 











-1 








1 


2 J 



\ 

-1 



0/ 



I)A,A-2 



dA,A-1 



5. Laplace operators 

The functions ([T]) are eigenfunctions of the Laplace operator A in variables x and y given relative to 
an orthogonal basis. The functions are also eigenfunctions of an algebraically independent degree four 
operator, say F, which is easily verified by direct computation. 

52 52 



A 



Qrf.2 Qy2 



r 



After the substitution of variables (16) and ( |17| ), the operators A and F are substantially transformed. 
When the operators are applied to the functions ([l]) expressed in the variables X, Y and U, V, W respec- 
tively, we obtain two operators for each family of polynomials ( 14 ), denoted here by A, F with following 



property: The polynomials corresponding to one of our four family of functions ( 14 ) are eigenfunctions 
of corresponding A,F [Ij. There are minute differences between them. Rather than explain the differ- 
ences, we write the two operators for all four cases. Note that the eigenvalues of A are zero for the 
lowest values of A and fi, admissible for each family of polynomials by the inequalities ([2]). 
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5.1. The operators A and F with eigenfunctions cosp^^^. 

Any polynomial P\^^{X, Y) of the cosp^^ family in variables ( 16 ) is the eigenfunction of the operator 

apa,m(^, y) = -(A' + ^i'')Px,M^ y) > (38) 

where 

~Kp^,^{x,y) = f(8-x2 + 2y)^ + (4x2-8y-2y2).^' 



'QX2 ^ 'QY^ 



Any polynomial P\^^{X, Y) of the cos^^ family in variables (16) is the eigenfunction of the fourth 



degree operator F = ^F, which is independent of the Laplace operator. 

fPx,^{X, Y) = X^fi^Px^X, Y) , (40) 



where 

TPxAX,Y) = (^(16 - + sy + y2)^ + (32X - + 16Xy + 2XY^)^^^ + 

(32y + i6y2 + 2y3 + lex^ + x^y^ - 4^^) g^2gy2 + (^2^^ - 8^^^ + i6xy2 + 2xy3)^^^+ 
(i6y2 - 4x2y2 + 8y3 + y^^— + (xy _ 4X)^ + (x^y - 24^^ + 4oy + 6y2 + 64) ^^^^^ + 
(-20x3 + 36xy + 7xy2 + 64X)— — + (64y - 2ox2y + 4oy2 + 6y3)— + y—^+ 

oXoY^ oY'^ oX'^ 

^^^^ " ^^^dm" + + + + + ^w) ^^^^ 

5.2. The operators A and F with eigenfunctions cos^^ /U. 

Any polynomial Px,fi{X,Y) of the cos^^^/C/ family in variables (16) is the eigenfunction of the 
operator A = (A + 1) 

A (Pa,m(^, Y)) = -(A2 + ^2 _ y) , (42) 

where 

A (P,,^(X, y)) = ((8 -X^ + 2Y)^ + (4X2 _ gy _ 2y2)^ 

+2(4X - XY)^^ - 3X A _ (4y + 8) A) P,,^(X, y) (43) 



Any polynomial Px^f^{X,Y) of the coSj^^^^/[/ family in variables (16) is the eigenfunction of the 
fourth degree operator F = :;pr^F, algebraically independent of the Laplace operator. 

TPx,,,{X, Y) = X^fi^PxjX, Y) , (44) 

where 



- / 8^ 

TPx,^,{X,Y) = f (16-4X2 + 8y + y2)— + (32X-8X3 + 16Xy + 2Xy2) 



+ 



dX^dY 



(32y + i6y2 + 2y3 + 16^^ + x^y^ _ 4^4) ___ + (32xy - 8x3y + i6xy2 + 2xy3)^^^^+ 



12 L. MOTLOCHOVA^ AND J. PATERA^ 



QY4. ^ ' QX'i ' ^ 

(-28x3 + 52xy + 9xy2 + 96X) ^^^^^ + (96y - 28x2y + 56y2 + sy^)^ + (2y + 4)^+ 



(5Xy - 8X)^^ + (-36X2 ^ ^ ^45^2 ^ gg)^ ^ ^ j p^^^(x,y) (45) 



5.3. The operators A and T with eigenfunctions sin^^^ /V. 

Any polynomial P\^^{X,Y) of the sin^^ family in variables (16) is the eigenfunction of the 

operator A = ^^jy (A + 2) 

A {Px,M, Y)) = -(A2 + ^?- 2)Px,^{X, Y) , (46) 

where 

A (Pa,m(X, y)) = f (8 - X2 + 2y)^ + (4X2 _ _ 2y2^ 5' 



9X2 'V ^QY^ 
_l Q Y _ 

9X9y 9x ay 



+2(4X - Xy) ^^^^^ - 3X— - 6y— ) Px,^,iX,Y) (47) 



Any polynomial Px^^{X, Y) of the sin^^ /!/ family in variables ( 16 ) is the eigenfunction of the fourth 
degree operator T = (F — 1), algebraically independent of the Laplace operator. 

TPxAX, Y) = (A2/i2 _ 1)P;,,^(X, y) , (48) 

where 



— 9^ d'^ 

TPx,^{X, Y) = ( (16 - 4X2 + gy + Y^)—- + (32X - SX^ + 16Xy + 2Xy2^ 



dX^ ' 'dX'^dY 

{32Y + 16y2 + 2y3 + 16X2 + X2y2 _ 4X4)— |^ ^ (32^y _ 8^3y ^ ^g^y2 ^ 2Xy3)— + 

aX^ar^ oXoY'^ 

(i6y2 - 4x2y2 + 8y3 + y^)— + (3xy - i2X)— ^ + (3x2y - 40x2 + 56y + ioy2 + 64)——+ 

a3 a3 

(-28x3 + 44xy + i3xy2 + 64X)— — ^ + (64y - 28x2y + 56y2 + ioy3^ 



dxdY^ ' 'ay3 

+(7y + x2 - 8)^ + (17Xy - 32X)^|^ + (-36X2 ^ 72y + 25y2 + 32)^ + 3X A 

+15y47)pxAX,Y) (49) 



9y, 

5.4. The operators A and T vi^ith eigenfunctions sin^^^ /W. 

Any polynomial Px^fj,{X,Y) of the sin^^^/PV^ family in variables (16) is the eigenfunction of the 
operator A = ^ (A + 5) 

A (Pa,m(^, Y)) = -(A2 + ^2 _ 5)p^^^(x, y) , (50) 

where 



„^v-, . - ^ vo - X2 + 2y)^ + (4X2 _ g^ _ 2Y')^ 



+2(4X - XY)^ - 5X A - (8y + 8)^ ) P,,,(X, Y) . (51) 



Any polynomial Px^fj,{X, Y) of the sin^.^ /y family in variables ( 16 ) is the eigenfunction of the degree 
four operator F = (F — 4), algebraically independent of the Laplace operator. 
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(52) 



where 



TPx,f,iX,Y) = i{l6-4:X^ + 8Y + Y'^)-^ + {32X-8X^ + l6XY + 2XY^) 



dX^dY 



+ 



(32y + i6y2 + 2y3 + lex^ + x^y^ _ 4^4)__^_ ^ (32j^y _ g^^ay ^ ;Lg^y2 ^ 2xy^)^^^^+ 

fl4 a3 ^3 

(16y2 _ 4x2y2 + 8y3 + y4)_ ^ (3^y _ 12X)-^ + (SX^y - 48^2 + 72y + 12y2 + 96) ^^^^^ 



(-36X3 ^ gQ^y ^ 15Xy2 + 96X)-— — + (96y - 36X2y + 72y2 + I2y3)— — 

'ax2 



+{X^ + lOy + 4)^ + (23Xy - 32X)— 1^ + (-64X2 ^ -^^Qy ^ 3gy2 + gg) + 5^ ^ 



ax 



+(32y + 32)^)PA,M(^,n 



(53) 



6. Continuous orthogonality of the four families of polynomials 

Continuous orthogonality of functions within each of the famihes cos^ , sin^ has been shown |141 [T5] . 
Recah that the functions of the same family are orthogonal when integrated over the triangle F (see Q). 



Our task here is to reformulate the orthogonality relations after the substitution of variable (16). The 



substitution transforms the domain of orthogonality integration F into its image ^ shown in Fig (6.1). 
Practically one needs to calculate the Jacobian of the change of integration variables. 



The Weyl formula |T3] for the integration of class functions over a compact simple Lie group G 
reduces the integration to the maximal torus T of G and independently provides information about the 
Jacobian. Up to some normalization conventions, it is the function 5^ of G. Our functions are W- 
invariant class functions, therefore integration is reduced to the fundamental region of the affine Weyl 



group W on T. These are commonly drawn conclusions for the functions (10) and (11). Corresponding 



conclusions for functions G and 5"+ of (12) and (13) are made only here. 



6.1. Orthogonality of the polynomials of . 

Consider any g\^cy^ Ax,y) and g\i^'cfy ,Ax,y) and their respective polynomials P^./x ^^'^ -^A',/. 



Computing directly the absolute value of the Jacobian 



D{X,Y) 
D(x,y) 



, we have 



J :-- 



DjXX) 
D(x,y) 



^7r2^4s(2,-^)(a;,y)S(2,i)(^'y) 



'(4,2) "-(4,0) 



'^(2,0) + 2 



^(3,3) + ^'^(3,1) 



^(1,1) 



27r^V |(-4X4 + 16X2 + 24X2y - 64y - 32y2 - AY^ + X2y2) 

7rV(^^ - 4y)((y + 4)2 - 4x2) = ^2 ^jjy^ ^ ^2 ^ 



(54) 



where U, V, and W are defined in (17). The fact that the Jacobian can be factorized in this way is of 
importance in the sequel. 



1 /-i 



Jy 





J-^PxA^^ Y)P^'y {X, Y) dX dY 



PxAX,Y)Py^^,{X,Y) 
4) V(^2_4y)((y + 4)2_4X2) 



dXdY = g\f,5xySi 



(55) 



14 



L. MOTLOCHOVAi AND J. PATERA^ 



V = ^^^^^ 



1 X 




4 X 



F = -2(X + 2 



Figure 1. Fundamental domain F in terms of the variables x,y, and its image after 



the change of variable (16) to X, Y . The three equations defining the boundary of 5^ are 
shown. 



6.2. Orthogonality of the polynomials of c /U. 

Consider any gx^^cV^ Ax, y) and gy^/cVy ,Jx, y) and their respective polynomials UPx^^ and UP\i^^ 



Multiplying c^^ and c^^, , one needs to know 



= 4c^i,0) • C(i,o) = 2cJ,o) + 4 - 4c+ = - 4y , 



which is the first of the two factors of the Jacobian (|54|). 
"1 ri 



Jy 



4 /"fi^+a /.4 /.4 



4J-iy-2 Jix2J-4^ 



I f /■|y+2 /.4 r4 



Pa,^(X, Y)Px',^,'{X, Y)U'J-' dX dY 



X2 - AY 



(56) 



6.3. Orthogonality of the polynomials of s'^ /V. 

Consider any and s^^, ^,^(x,y) and their respective polynomials VPx^j^ and VPX'fi'. 

Multiplying Sp^ and Sp^, one needs to know 

= 4^1,1) • <i,i) = 2c+ 2) - 4cJ,o) + 4 = (y + 4)2 - 4X2 ^ 



which is the second of the two factors of the Jacobian (|54|) 

»1 rl 

9^i^st,f.)9\'f.'S^y^^,^dxdy = 



Jy 



4 /"l^+S ^4 />4 



4J-|y-2 Jljs:2j-4^ 



Pa,m(^, ^) A'y Y)V^J-^ dX dY 



1 



4 ri ^4 



vr" \ .1-4 J-\Y-2 J|X2 J-4, 



(y + 4)2 - 4X2 



(57) 



6.4. Orthogonality of the polynomials of the antisymmetric sine functions. 

Consider any s^^^(x,y) and s^, ^,^(x,y) and their respective polynomials WPx/j, and WPy^'- 

Multiplying and s^, one needs to know 

^' = 4S(2,l)-^(2,l) = ^'^'' 



which is up to the constant vr ^ the square of the Jacobian (54). 
rl rl 



'0 Jy 



^A,^S(A,^)5AV'S(A'y)'^^'^^ 



4 r^Y+2 ri rl 



Y~2 Jlx'^J-A 



Px,^,{XX)Px',^.'{XX)W^J-UXdY 
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= 9\^lh\'5^J,^' (58) 
7. Identities of generalized trigonometric functions 

Identities verified by the common trigonometric functions of one variable are numerous, well known, 
and often useful. Many more relations can be found for the two variable generalizations ([T]) of the 
trigonometric functions. In this section, there are written relations for general values of parameters and 
the only condition on the parameters is that they are integers and distinct, i.e. they need not to satisfy 
the inequalities ([2]). Relations, obtained when some of the parameters coincide, are simpler. Numerous 
such relations are found in Appendix 1. Analogs of known identities of one variable trigonometric 
functions are more easily recognized there. 

Relations are arranged according to the type of functions ([T]) appearing on the right hand side of 
each equality. 

7.1. Generic relations for c"*". 

^A,m)^A',m') + ^(A,M)^(A',/.') + ^(A,/.)S(A',M') + ^(A,m)^(A',m') = ^(A+A',M+M') + ^A-A',^*-^') ^^^'> 

^A,m)4a',m') + ^(A,M)^(A',/.') - ^(A,^)«(A',At') - ■\\,^,f{\',^,') = C(A+A',M-M') + ^A-A'.^^+m') (^°) 

'=(A,m)'^(A',m') ~ '^(A,/.)C{A',/.') + ■5(A,^)^(A',M') ~ ■\x,^,f{\',^l') = ^(A+M'.M+A') + C(A-^',A.-A') (^1) 

'^tx.u)'^t\'.u') ~ '^(X.a)^(X'.u') ~ ^tx.u)^tx'.u') + ^ ( X.u)^ ( X' .u') = ^tx4- u' . u- X' ) + ^tx-u' .u+X') i^"^) 



7.2. Generic relations for c . 

2 {^t^,fif(X',fi') + ^(A,M)'^(A'y)) ^ '^{A+A',/i+At') + '^(A-A',M-/i') ~ '^{A+m'.m+A') ~ ^{X-^i' ,fi-X') (^'^) 
2 ( '^tx.ll^'^(X'.ll'^ ~ ^tx.ii.)^(X'ii') ) ~ '^(X+X'.a-a') + '^fA-A'.H+uM ~ ^(X+ii'a-X') ~ X-a'.ii+X') (^^) 



''(X,^l)'-{X' ■'(A,m)'^(A',aj'); ~ {A+A',/i-At') (A-A',M+/i') ^{X+^l',^^-X') "-(A-At'./i+A') 

7.3. Generic relations for s~^. 

2 (*(A,m)'^(A',m') + ^{X,fif{X',fi')) = ^(A+A',M+A«') + ^(A+A',M-/^') + ^(A-A'.aa+a*') + ^(A-A',^*-/.') (^5) 
^ ('^(A,m)'^(A',m') ~ '^(A,/x)^{A',/x')) " ^(A+Af ,A»+A') + '^(A+m',/^-A') + ^(A-/.',a»+A') + '^(A-/.',/.-A') (^6) 

7.4. Generic relations for s . 

^ ('^(A,Ai)'^(A',At') + ^iX,fif{X'y)) = ^(A+A',M+A«') + '^(A+A',^-/^') + ^ {X-X' + ■^(A-A',^*-/.') (^'^) 
2 ('^(A,m)'^(A',m') ~ '5(A,/.)^(A',/x')) " ^(A+M',|»+A') + ■5(A+/x',/.-A') + ^(A-^'./^+A') + «(A-/.',/.-A') (^^) 

8. Conclusions 

The explicit relation between the functions studied here as 'generalizations of 2-dimensional trigono- 
metric functions' ([T]), and the same functions after the linear substitution of variables ([9]), called the 
orbit functions of the simple Lie algebra or group of type C2 (see (10)-(13)), has important practical 
and theoretical implications. Indeed, the functions are simpler in their 'trigonometric' form ([T]) than in 
their C2 form, and are given in variables that refer to an orthogonal basis in M?. The Lie theoretical 
perspective on the functions allows one to discretize the orbit functions and carry over the discretization 
to the four families of C2 polynomials. 

It is easy to recognize that among the four families (14) the last one is formed by the character of 
finite dimensional irreducible representations of C2. Perhaps even more significant are the implications 
of the fact that the second and third families of functions in ( 14 ) have not been encountered so far in 
the literature or in the Lie theory in general. 

As a result of this paper, the number of families of orthogonal polynomials based on C2 that can be 
built has grown from three (C-, S-, and i?- functions of [23]) to four (C^ and S^), and to another six 
ii^-like combinations + . 
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Generalization of these facts to all simple Lie groups with two lengths of simple roots is imminent. 

Interesting and undoubtedly useful are the new functions C~ and due to their 'hybrid' behaviour 
at the boundary of the region F, i.e. they are symmetric on some sides of the simplex F, skew-symmetric 
on the other sides. 

There is another related but different point of view on the polynomials we constructed in this paper. 



Tables in Appendix 2 show examples of the polynomials in the form (14), whose variables X and 



Y are the two lowest C-functions (16). Alternatively, one may choose as the polynomial variables the 



characters Xui and Xui2 of the C2 representations of the lowest dimension, namely 4 and 5. The character 



variables are then related to X and y of (16) as follows, 



Xi^l=X , Xa;2 = ^ + 1 • 

The characters have all of the properties required to build the polynomials, including continuous and 
discrete orthogonality [23]. Often the characters are avoided in applications because, unlike orbit 
functions, they grow without limits for higher representations. However, on the level of polynomial 
exploitation, that difference is not visible. Hence the alternative can be rather useful |10j . 
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l. motlochovai and j. patera^ 
Appendix 1 



In addition to the generic recurrence relations of Sec. [7| here we present simphfied relations obtained 
when some of the parametres of (59)-(68) coincide. 



Relations for 



A' = 


A 
















''(Ay) 




+ *{A,m)^(A,/.') 


+ ^(A,m)'^{A,/.') 


= '^(2A,/i+M') '^(0,M- 








''(A.A.') 


+ C(A,m)C(A,m') 


*(A,m)^(A,/.') 




= '^(2A,/i-M') + '^(0,M^ 








^(Ay) 


- C(A,m)^(A,m') 


+ *(A,m)^(A,/.') 


~ ^{X,^^f(X,t,') 


= '^(A+/i',M+A) + ^(A- 


-m'.m- 


-A) 








*{A,m)^(A,/.') 


+ ^(A,m)'^{A,/.') 


= '^(A+/i',At-A) + ^{A- 




^A) 



A' 

4. 

c 
c 
c 
c 



(a,m: 

+ 

(a,m: 
+ 

(a,m: 



(a,m: 



A 



+ 

(A,A 
+ 

(A,A 



'(A,A 
-(A.A 



'(A,A 
'(A,A 



+ 

(A,A 
+ 

(A,A 



^(A,/.) +C(A,m)^(A,m) 
C^A,/.) +C{A,m)^(A,m) 



'(A,m)-'{A,m) 



'(A,m)^(A,m) 



^(A,m)'^{A,/.) ■5{A,/.)'5(A,m) 



^(2A,2/^) + ^ 
'^{2A,0) + '^(0.2^) 



•^(A./i) ^{x,^lf{x,^Ji) + *(a,m)'^{a,/x) '^{a,/x)^(a,m) ~ ^(A+^.A+^t) + C(^_a,^_a) 



1A,/.) 



^{A,/.)C(A,m) 



^{x,^Jif{x,^l) + ^{x,^lf(x,^Ji) 



2c+ 

(A+^i,A-/i) 



'(A',Ai') ^ ''(A,A)''{A',Ai') - ^(A+A',A+Ai') ^ "-(A-A'.A-^t') 



,n + c; 



{X',^l') '^{X,Xr{X',/i') ~ ^{X+X',X-/i') ^ "-(A-A'.A+^t') 

A, // = A' 

'^(A',A') + '^(A,A)'^(A',A') ~ '^(A+A',A+A') + "^(A-A'.A-A') 



'^{A',A') ^(A,A)'^(A',A') 



A' = A 



2c 



(A+A',A-A') 



^{A,m') + '^{A,A)^(A,m') ~ ^(2X,X+/i') + "^(O.A-mO 
^{A,m') ~ '^{A,A)^(A,m') ^ '^(2A,A-/x') + "^(O.A+m') 



/U = A 



'^{A,A) + '^(A,A)'^(A,A) 
'^(A,A) ~ '^(A,A)'^(A,A) 



'^(2A,2A) + ^ 



2c: 



(2A,0) 



Relations for c 



A' = A 



2|C(A,^)C(A,^') +^(A,m)*(A,/.') 
2 'c^A,Ai)C(A,M') - ^(A,m)^(A,/.') 



^{2X,^l+fl') ^{A+M',At+A) '^(A-At',M-A) 



'^(2A,fi-/i') '^{A+At',At-A) '^(A-At',^+A) 
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A' = A, fJ'' = IJ- 

2 (cfA,M)^(A,M) + *(A,m)*(A,m)) = ^(2A,2m) 

2 (^t^,t,fi\n) - ^(A,,.)^(A,m)) = ^(2A,0) + ^(0,2/.) 

/X = A 

'^{A,A)'^(A',/i') + ^t^,\)^{X',fi') ^ "^(A+A'.A+^t') + '^(A-A',A-^') 

^t>^,xf{x'y) ~ '*(a,a)'^(a',m') " '^(a+a',a-m') + '^(a-a',a+,.') 
/7, = A' = A 

'^(A,A)'^(A,^') + ^t^,xf(X,tM') ^ '^(2A,A+^j') 



^(A,A)'^(A,/.') ^{X,xfiX,n') - ^(2A,A-m') 



Relations for s+. 

V = A 



2 (^^(A,/.)'^(A,/.') + *(A,^^)^(A,^')j ~ ^i2X,^l+^^'} + '^(2A,m-m') 

2 (■S(A,m)'^(A,m') ~ ■*(A,m)'^(A,/.')) = '*(A+M',M+A) + ■*(A+/.',m-A) + •S(A-/.',/.+A) + •S(A-,x',,.-A) 

A' = A, IJ'' = 

2 («(A,m)^(A,m) + *(A,m)^(A,m)) = *5a,2m) 

^ (*{A,m)'^(A,//) *(A,/i)'^(A,;ti)) ~ *(A+/i,A+;i) *(A-//,A-ai) 

/Lt = A 

2*(A,A)'^(A',m') ^ *(A+A',A+/i') + ^(A+A',A-;i') + ^(A-A',A+/i') + *(A-A',A-|u') 

fj, = X, = X' 

^*(A,A)'^(A',A') ^ '^(A+A',A+A') + ^*(A+A',A-A') + ■^(A-A'.A-A') 

n = x' = x 

^^t>^,xft>^,n') " ^5a,A+;.') + *(2A,A-m') 

^ = yti' = a' = A 

2«(A,A)C(A,A) = *(2A,2m) 

Relations for s~. 

X' = X 

2 (*(A,m)'^(A,m') + *(A,m)'^(A,m')) = *(2A,M+M') + *(2A,/.-M') 

2 (*(A,m)'^(A,m') ~ *(A,m)'^(A,m')) " *(A+m',m+A) + «(A+m',m-A) + *(A-m',M+A) + *(A-;u',M-A) 
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A' = A, 



/i = ^^ 



(A— /^,A+/i) 



/X = A 

^'^(A,A)'^(A',^i') ~ •^(A+A',A+/x') + ^{X+X',X~ii') + *(A-A',A+m') ~'~ "^(A-A',A-/i') 



2s 



A' = A 



(A,A)^(A,/.') 



^{2X,X+fj.') + ^(2A,A-/x') 



Appendix 2 



There are two tables shown for each of the four function of (14), one for each congruence class ([s]). 
The first table contains the functions of congruence class ^ = and the second table contains the 
functions with # = 1. 



c+ 


1 


Y 




y2 




y3 


X^ 


X^Y^ 


y4 


X^Y 




Y" 


■^6 






y-6 


X^Y 


X^Y^ 


X^Y" 


Y^ 


i(0,0) 


1 








































2(1,1) 





1 






































2(2,0) 


-4 


-2 


1 




































2(2,2) 


4 


4 


-2 


1 


































4(3,1) 





-6 





-2 


1 
































2(3,3) 





9 





6 


-3 


1 






























2(4,0) 


4 


8 


-4 


2 


-4 





1 




























4(4,2) 


-8 


-20 


10 


-12 


8 


-2 


-2 


1 


























2(4,4) 


4 


16 


-8 


20 


-8 


8 


2 


-4 


1 
























4(5,1) 





10 





10 


-5 


2 





-4 





1 






















4(5,3) 





-30 





-40 


20 


-16 





12 


-2 


-3 


1 




















2(5,5) 





25 





50 


-25 


35 





-20 


10 


5 


-5 


1 


















2(6,0) 


-4 


-18 


9 


-12 


24 


-2 


-6 


9 





-6 








1 
















4(6,2) 


8 


40 


-20 


42 


-48 


16 


12 


-24 


2 


12 


-4 





-2 


1 














4(6,4) 


-8 


-52 


26 


-100 


56 


-70 


-14 


51 


-20 


-12 


16 


-2 


2 


-4 


1 












2(6,6) 


4 


36 


-18 


105 


-48 


112 


12 


-72 


54 


12 


-36 


12 


-2 


9 


-6 


1 










4(7,1) 





-14 





-28 


14 


-14 





28 


-2 


-7 


9 








-6 








1 








4(7,3) 





42 





98 


-49 


70 





-84 


20 


21 


-34 


2 





18 


-4 





-3 


1 






4(7,5) 





-70 





-210 


105 


-224 





168 


-108 


-42 


94 


-24 





-30 


20 


-2 


5 


-5 


1 




2(7, 7) 





49 





196 


-98 


294 





-196 


210 


49 


-161 


77 





42 


-56 


14 


-7 


14 


-7 


1 



Table 2. 



The polynomials c^^ of the congruence class 7^ 



of degrees < 7. In the 



first column the labels (A,/i) are shown. The constants ^, 2 and 4 should be understood 
as multiplying the functions ^-^ , making all the entries of the table integer. 



c+ 


X 


XY 




XY'^ 




XY'-^ 


X" 




XY^ 


X^'Y 




XY^' 










2(1,0) 


1 
































4(2,1) 


-2 


1 






























2(3,0) 


-3 


-3 


1 




























4(3,2) 


6 


3 


-2 


1 


























4(4,1) 


2 


-4 





-3 


1 
























4(4,3) 


-6 


6 


2 


5 


-3 


1 






















2(5,0) 


5 


15 


-5 


5 


-5 





1 




















4(5,2) 


-10 


-25 


10 


-15 


10 


-3 


-2 


1 


















4(5,4) 


10 


10 


-10 


15 


-5 


7 


2 


-4 


1 
















4(6,1) 


-2 


9 





18 


-6 


5 





-5 





1 














4(6,3) 


6 


—21 


—2 


—45 


18 


—21 





15 


—3 


—3 


1 












4(6,5) 


-10 


15 


10 


35 


-20 


28 


-2 


-16 


9 


5 


-5 


1 










2(7,0) 


-7 


-42 


14 


-35 


35 


-7 


-7 


14 





-7 








1 








4(7,2) 


14 


77 


-28 


84 


-70 


35 


14 


-35 


5 


14 


-5 





-2 


1 






4(7,4) 


-14 


-56 


28 


-105 


63 


-84 


-14 


56 


-27 


-14 


20 


-3 


2 


-4 


1 




4(7,6) 


14 


21 


-28 


70 


-28 


84 


14 


-56 


45 


7 


-31 


11 


-2 


9 


-6 


1 



Table 3. The polynomials of the congruence class # = 1 of degrees < 7. In the 

first column the labels (A, fi) are shown. The constants 2 and 4 should be understood as 
multiplying the functions ^ , making all the entries of the table integer. 





X 


XY 


X'' 


XY'^ 


X-^Y 




X" 




XY^ 


X-'Y 


^3y3 


XY'' 


2(2,0) 


1 
























4(3,1) 





1 






















2(4,0) 


-4 


-2 


1 




















4(4,2) 


6 


4 


-2 


1 


















4(5,1) 





-5 





-2 


1 
















4(5,3) 





10 





6 


-3 


1 














2(6,0) 


9 


12 


-6 


3 


-4 





1 












4(6,2) 


-16 


-24 


12 


-12 


8 


-2 


-2 


1 










4(6,4) 


10 


20 


-10 


21 


-8 


8 


2 


-4 


1 








4(7,1) 





14 





14 


-7 


3 





-4 





1 






4(7,3) 





-35 





-42 


21 


-16 





12 


-2 


-3 


1 




4(7,5) 





35 





56 


-28 


36 





-20 


10 


5 


-5 


1 



Table 4. The polynomials i7 ^c^^^^ of the congruence class # = of degrees < 6. In 
the first column the labels (A, ji) are shown. The constants 2 and 4 should be understood 
as multiplying the functions U~^c7. x, making all the entries of the table integer. 
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■ 


1 


Y 




y2 




y3 


X^ 




y4 


X'^Y 














2(1,0) 


1 
































4(2,1) 


2 


1 






























2(3,0) 


-3 


-1 


1 




























4(3,2) 


6 


5 


-2 


1 


























4(4,1) 


-2 


-4 





-1 


1 
























4(4,3) 


6 


14 


-2 


7 


-3 


1 






















2(5,0) 


5 


5 


-5 


1 


-3 





1 




















4(5,2) 


-10 


-15 


10 


-7 


6 


-1 


-2 


1 


















4(5,4) 


10 


30 


-10 


27 


-11 


9 


2 


-4 


1 
















4(6,1) 


2 


9 





6 


-6 


1 





-3 





1 














4(6,3) 


—6 


—29 


2 


—27 


18 


—9 





9 


— 1 


—3 


1 












4(6,5) 


10 


55 


-10 


77 


-36 


44 


2 


-24 


11 


5 


-5 


1 










2(7.0) 


-7 


-14 


14 


-7 


21 


-1 


-7 


6 





-5 








1 








4(7,2) 


14 


35 


-28 


28 


-42 


9 


14 


-19 


1 


10 


-3 





-2 


1 






4(7,4) 


-14 


-56 


28 


-77 


49 


-44 


-14 


40 


-11 


-10 


12 


-1 


2 


-4 


1 




4(7,6) 


14 


91 


-28 


182 


-84 


156 


14 


-96 


65 


17 


-41 


13 


-2 


9 


-6 


1 



Table 5. Polynomials U C(^Xfj,) congruence class # = 1 of degrees < 6. La- 

bels (A, n) are given in the first column. Constants 2 and 4 should be understood as 
multiplying functions U~^c7. ^ making all entries of the table integer. 





1 


Y 


X-' 


y2 


X'^Y 


y3 


X^ 


^2y2 


Y* 




^2y3 


y5 




^4y2 


X^Y* 


y6 


2(1,1) 


1 
































2(2,2) 





1 






























4(3,1) 


-2 


-2 


1 




























2(3,3) 


1 


2 


-1 


1 


























4(4, 2) 





-4 





-2 


1 
























2(4,4) 





4 





4 


-2 


1 






















4(5,1) 


2 


6 


-3 


2 


-4 





1 




















4(5,3) 


-2 


-8 


4 


-8 


4 


-2 


-1 


1 


















2(5,5) 


1 


6 


-3 


11 


-4 


6 


1 


-3 


1 
















4(6,2) 





6 





8 


-4 


2 





-4 





1 














4(6,4) 





-12 





-22 


11 


-12 





8 


-2 


-2 


1 












2(6,6) 





9 





24 


-12 


22 





-12 


8 


3 


-4 


1 










4(7,1) 


-2 


-12 


6 


-10 


20 


-2 


-5 


9 





-6 








1 








4(7, 3) 


2 


14 


-7 


22 


-20 


12 


5 


-14 


2 


6 


-4 





-1 


1 






4(7, 5) 


-2 


-18 


9 


-48 


24 


-44 


-6 


30 


-16 


-6 


12 


-2 


1 


-3 


1 




2(7, 7) 


1 


12 


-6 


46 


-20 


62 


5 


-39 


37 


6 


-24 


10 


-1 


6 


-5 


1 



Table 6. Ppolynomials V ^s'l^Xfx) congruence class # = of degrees < 6. La- 

bels (A, //) are given in the first column. Constants 2 and 4 should be understood as 
multiplying the functions x, making all entries of the table integer. 



4- 


X 


XY 






X'-^Y 


XY'-'' 


X" 


X'''Y'^ 










4(2,1) 


1 
























4(3,2) 


-1 


1 






















4(4, 1) 


-2 


-3 


1 




















4(4, 3) 


2 


1 


-1 


1 


















4(5,2) 


1 


-3 





-3 


1 
















4(5,4) 


-2 


3 


1 


3 


-2 


1 














4(6,1) 


3 


12 


—4 


5 


-5 





1 












4(6,3) 


-3 


-9 


4 


-9 


5 


-3 


-1 


1 










4(6,5) 


3 


3 


-4 


8 


-2 


5 


1 


-3 


1 








4(7,2) 


-1 


6 





15 


-5 


5 





-5 





1 






4(7,4) 


2 


-9 


-1 


-24 


10 


-15 





10 


-3 


-2 


1 




4(7,6) 


-3 


6 


4 


16 


-10 


17 


-1 


-9 


7 


3 


-4 


1 



Table 7. Polynomials V s'^xix) congruence class # = 1 of degrees < 6. Labels 

(A, jj) are given in the first column. The constant 4 should be understood as multiplying 
the functions V~^s1'. ^, making all entries of the table integer. 





X 


XY 


X^ 


xy2 


X''Y 


xy-^ 






xy* 


4(3, 1) 


1 


















4(4,2) 





1 
















4(5, 1) 


-3 


-2 


1 














4(5,3) 


2 


2 


-1 


1 












4(6,2) 





-4 





-2 


1 










4(6,4) 





5 





4 


-2 


1 








4(7,1) 


6 


10 


-5 


3 


-4 





1 






4(7,3) 


-5 


-10 


5 


-8 


4 


-2 


-1 


1 




4(7,5) 


3 


8 


-4 


12 


-4 


6 


1 


-3 


1 



Table 8. Polynomials W s^^^^ of the congruence class # = of degrees < 5. Labels 
(A, jj) are given in the first column. The constant 4 should be understood as multiplying 
the functions W~^s7. x, making all the entries of the table integer. 



w" 


1 


y 


X' 


y2 


x^y 


Y' 


X^ 




y4 


x^y 




y5 


4(2,1) 


1 
























4(3,2) 


1 


1 






















4(4,1) 


-2 


-1 


1 




















4(4,3) 


2 


3 


-1 


1 


















4(5,2) 


-1 


-3 





-1 


1 
















4(5,4) 


2 


7 


-1 


5 


-2 


1 














4(6,1) 


3 


4 


-4 


1 


-3 





1 












4(6,3) 


-3 


-7 


4 


-5 


3 


-1 


-1 


1 










4(6,5) 


3 


13 


-4 


16 


-6 


7 


1 


-3 


1 








4(7,2) 


1 


6 





5 


-5 


1 





-3 





1 






4(7,4) 


-2 


-13 


1 


-16 


10 


-7 





6 


-1 


-2 


1 




4(7,6) 


3 


22 


-4 


40 


-18 


29 


1 


-15 


9 


3 


-4 


1 



Table 9. Polynomials W s^^^^^oi the congruence class 7^ = 1 of degrees < 5. Labels 
(A, /x) are given in the first column. The constant 4 should be understood as multiplying 
the functions W~^s7y^ making all the entries of the table integer. 



